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Abstract 

In this paper we investigate the parameter estimation of the fiber lay-down process 
in the production of nonwovens. The parameter estimation is based on the mass per unit 
area data, which is available at least on an industrial scale. We introduce a stochastic 
model to represent the fiber lay-down and through the model's parameters we characterize 
this fiber lay-down. Based on the occupation time, which is the equivalent quantity for 
the mass per unit area in the context of stochastic dynamical systems, an optimization 
procedure is formulated that estimates the parameters of the model. The optimization 
procedure is tested using occupation time data given by Monte-Garlo simulations. The 
feasibility of the optimization procedure on an industrial level is tested using the fiber 
paths simulated by the industrial software FYDIST. 

1 Introduction 

Nonwoven materials or fleece are webs of long flexible fibers that are used for composite mate- 
rials, e.g. filters, as well as in the hygiene and textile industries. They are produced in melt- 
spinning operations: hundreds of individual, endless fibers are obtained by the continuous 
extrusion of a molten polymer through narrow nozzles, which are densely and equidistantly 
placed in a row at a spinning beam. The viscous or viscoelastic fibers are stretched and spun 
until they solidify due to cooling air streams. Before the elastic fibers lay down on a moving 
conveyor belt to form a web, they become entangled and form loops due to highly turbulent 
air flows. Figure [T] shows at a microscopic level the webs formed on the conveyor belt. The 
homogeneity and load capacity of the fiber web are the most important textile properties for 
quality assessment of industrial nonwoven fabrics. The optimization and control of the fieece 
quality require modeling and simulation of the fiber dynamics and the lay-down. 
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Figure 1: Forms of fiber webs (nonwovens). 



There are two classes of approaches to model the fiber lay-down process. The first class 
uses microscopic details to model the lay-down of the fibers. Here, each fiber is seen as an 
elastic beam. The software FYDIST, developed by the Fraunhofer ITWM, Kaiserslautern, 
Germany, uses such models to describe the fiber lay-down. Since the motion of each fiber 
is simulated using the physics on a microscopic level, the behavior of the simulated fibers is 
quite close to the real industrial fibers. Nevertheless, due to the large number of microscopic 
details included, the numerical computations are highly time consuming. 

The second approach is based on macroscopic, quantitative description. Here, the lay- 
down is modeled stochastically, i.e. the models are consist of stochastic differential equations 
with a certain set of parameters. Since this quantitative approach does not use fine details of 
the fiber lay-down, it allows for fast numerical simulations. 

Available data to judge the quality, at least on the industrial scale, are usually the mass per 
unit area of the fleece. Figure [2] (right) shows the mass distribution per area of a fieece made 
by a single fiber (left). Since the mass distribution is an averaged quantity, the microscopic 
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Figure 2: A fiber path (left) laid on a conveyor belt and the mass distribution of the fiber 
(right). 

details of the fibers do not play a significant role. Nevertheless, the quantitative, macroscopic 
approach provides sufficient details to optimize fibers based on mass distributions. Moreover, 
due to fast simulations the stochastic models increase the efficiency of the optimization. Thus, 
in order to identify the parameter of the fiber lay-down we use the macroscopic approach. 
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A stochastic model for the fiber deposition in the nonwoven production was proposed and 
analyzed in Ref. [21 [3]. In particular, the stochastic model proposed in Ref. [2J, which repre- 
sents the fiber deposition on a moving belt, is constituted by nonlinear stochastic differential 
equations (SDE). 

The hydrodynamic scaling limit of the resulting stochastic process of the model is given by 
an Ornstein-Uhlenbeck process with moving mean. The aim of this paper is to determine the 
parameters of the Ornstein-Uhlenbeck process with moving mean from available mass per 
unit area data, i.e. the occupation time. 

The paper is organized as follows: In Section [2] we introduce the fiber lay-down process as 
a simplified model for the fiber deposition. Furthermore, we define the occupation time and 
present an analytical expression of the expected occupation time. An optimization procedure 
to estimate the model parameters from available occupation times is presented along with 
numerical experiments in Section [3} Finally, we draw some conclusions and give an outlook 
to open questions. 



2 Model and theory 

Summarizing [2J , we model the fiber lay-down process on a moving conveyor belt by a stochas- 
tic process Y := {Y^^ ,Y^)f~^Q E satisfying the stochastic differential equations 

dY} = X(6i + Kt- y/) dt + (Ji dW}, 

(1) 

dy/ = \{02 - Yt^) dt + (72 dW^, 

with initial condition Yq = {Yf^,Y^)^ G M^. The stiffness A > governs the deterministic part 
of the fiber deposition, while a standard two-dimensional Brownian motion W = {W^, W^)fyQ 
and the diffusion parameters di, (T2 G M govern the stochastic part. The parameter k > is 
the belt speed and (6*1, 6*2)^ G is the reference position of the fiber lay-down. 
The fiber lay-down process is always centered to its moving position (^i + Kt, 92). Assuming 
the deterministic forces are symmetric around the reference position of the fiber lay-down, 
we use the same stiffness coefficient A in both equations, i.e. an isotropic fiber lay-down. For 
sake of simplicity, we set 9i = 62 = 0. Then, the fiber lay-down model reads as 

dY^ = -\{Y^ - Kt) dt + ai dWl, 

dY^ = -\Y^dt + a2 dW^ . 

The random variable Yt models the deposition point of an individual fiber on the fleece. If 
we follow the random variable over a time interval [0, T] for T > 0, we obtain the path of an 
individual fiber. The fiber lay-down process Y on a moving belt is characterized by the drift 
and diffusion parameters A, ai and o"2, together with the belt speed k. Assuming the belt 
speed K as a known information, the process Y given by the parameters A, ai and ct2 may be 
denoted as yA,cri,(T2- 

To introduce the mathematical analogue of the mass per unit area we need the following 
definition. 

Definition 2.1 (Occupation time). Let T > and consider a rectangle V := [ai,6i] x 
[02)^2] C M^. The occupation time Mx>,t is defined as 

Mv,T{y) ■= riv{Yt)dt. 
Jo 
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Here, lx> denotes the indicator function of the rectangle T). 



Remark 2.1. The occupation time is a random variable itself. It models the time, the random 
process spends inside the rectangle T> during the time interval [0,T]. In terms of our physical 
model for the nonwoven production, the occupation time can be interpreted as the mass of 
fiber material deposited inside D, i.e. the mass per area of the final fleece. This quantity is 
easily accessible even on the scale of industrial production and hence it will serve as the input 
to our parameter estimation problem. 

The following theorem can be proven using the techniques of white noise analysis, see [1] 
for the one-dimensional case. 

Theorem 2.1 (Occupation time). Xei Mjq^ ;,^]^ [02,62], t(^) the occupation time of the fiber 
lay-down process Y given by ^ . Then, its expectation is given by 

1 }f^^J VXib^-f{Xt-l + e-'^)) \_^^J VX{a,- f{Xt-l + e-'^))\\ 




0-2^1 — «~2At 




X erf -erf (3) 



In the following we use E\ax^hx\Y.{a2,h2\ k, tI-^j cti, (72) to denote the expected occupation time 
of the process 1a,cti,o-2- 



3 Numerics 



3.1 Monte— Carlo computation of the expected occupation time 

A Monte-Carlo method that approximates the expected occupation time of the fiber lay-down 
process Y consists of the following steps 

• approximate numerically the process Y given by the SDE ([2]), 

• approximate the expected occupation time by considering a sufficiently large number of 
sample paths of a two-dimensional Brownian motion W . 

To approximate the process numerically we use the standard Euler-Maruyama method up to 
time horizon [0, T], see [5]- Let us denote the Euler-Maruyama approximation on the time 
grid {0 = To < Ti < • • • < Tfc = r} by Y^^^ . Figure [s] shows two sample paths of the process 
Y . Both sample paths are plotted with respect to the Cartesian coordinate system defined 
by the nozzle position and the moving direction of the conveyor belt. The axis that passes 
through the nozzle position and parallel to the moving direction of the conveyor belt is called 
the central axis of the fiber lay-down. 

Let i denote the mapping implicitly defined by the stochastic differential equations ([2]), 
i.e. Y = ^{W) for a given sample path of the Brownian motion W. Let ^^^^ be the mapping 
implicitly given by the Euler-Maruyama scheme, i.e. Y^'^^ = ^^'^^W). Thus, computing the 
integral 

M[a,,b,Ma2M,-,Ti^'^'^(W))= [ lK,Hx[a2.62](^'^'HW^))di, 

Jo 
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Figure 3: Sample path of the process Y given by the parameters A = 1, cii = 1 and cr2 = 1. 
The time horizon and the belt speed are T = 30 and k = 0.5, respectively. 

we get the occupation time of Y^'^^ The Monte-Carlo approximation for the expected occu- 
pation time EiMia^^biMaiMl.i^.Tiy)) is given by 

1 ^ 

5^(I>«(H^)) = j^Y.^M]>^MU,t{^^'Hw,)), (4) 
i=i 

where Wi , W2 , • • . , Wn are independent sample paths of the two-dimensional Brownian motion 
W. Due to the law of large numbers we know that Sn{Y^^'') — )• E{M^ni,bi]xla2,b2], k,t{Y)) ^ 
N 00. For more details see [6j and references therein. 

Table [T] lists computations of the expected occupation times for the process Y corre- 
sponding to the area [3, 15] x [—1,1] for different parameters A, ai, 02-, belt speeds k and time 
horizons [0, T]. The Monte-Carlo method uses 10000 independent sample paths. 



Parameters 

Time T 


K = 1 
A, CTl, (72 = 1 


K = 2 
A, 0"1, (T2 = 1 


K = 1 
A, (Tl, (72 = 2 


K = 2 
A, (Tl, CJ2 = 2 


T = 7 by MC 


2.5465 


3.8793 


2.3899 


3.4246 


T = 7 by ( 


3 




2.5484 


3.8742 


2.3910 


3.4207 


T = 30 by MC 


10.1382 


5.1327 


8.1910 


4.1021 


T = 30 by 


3 


) 


10.1327 


5.1380 


8.1938 


4.1063 


T = 50 by MC 


10.1382 


5.1327 


8.1910 


4.1021 


T = 50 by 


3 


) 


10.1327 


5.1380 


8.1938 


4.1063 



Table 1: Expected occupation times of Y computed by the Monte-Carlo (MC) method and 
the formula ([s]) . 

Note that for given parameters A, ai and (T2 the expected occupation times corresponding 
to the time horizon T = 30 do not differ from those corresponding to the time horizon T = 50. 
This observation is due to the movement of the belt in the lay-down process. Note that in 
cases T = 30 and T = 50 the Monte-Carlo method uses the fiber sample path corresponding 
to random sequences generated by a fixed set of seeds. 

Figure [4] compares the computation of expected occupation times computed by the ana- 
lytical formula (Is]) and the Monte-Carlo method with 10000 sample paths. 
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Figure 4: Comparison of the expected occupation times -Ec, r(A, ui, (J2) computed by the 
formula ([s]) and the Monte-Carlo method for the domain T) = [1, 3] x [—1.5, 1 — 5], belt speed 
HI = 2 and time horizon T = 20. 



To analyze how the expected occupation time of a given process Y is distributed around 
the central axis of the fiber lay-down, we proceed as follows: 

Let us consider a fiber lay-down process Y given by the parameters A = 1.5, o"! = 1 and 
(T2 = 1. We fix the belt speed k = 1, time horizon T = 20 and consider an area given by 
[5, 15] X [—1.8, 1.8], which is symmetric around the central axis. We subdivide this area into 
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Figure 5: The grid that divides the area [5, 15] x [—1.8, 1.8] (left) and the expected occupation 
time distribution (right) with respected to the cell areas given by the grid. 



a grid of 25 x 15 rectangular sub areas (cells) of equal sized, see Figure [s] (left). Next, we 
compute the expected occupation times corresponding to each cell using the formula ([s]) . On 
the right hand side of Figure [5] we plot the corresponding expected occupation times. Note 
that the occupation time decreases as we move away from the central axis y2 = of the lay- 
down process. Increasing the parameter a2 or decreasing the parameter A leads to a flater 
and broader distribution of the expected occupation time with respect to the y^-coordinate. 
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3.2 An optimization procedure for the parameter estimation 



As we have aheady mentioned, the fiber lay-down process Y is characterized by the parame- 
ters A, fJi, (T2 and the belt speed k. Using this relation inversely, we estimate the parameters 
A, cJi and cj2 based on a given set of expected occupation times. Including the belt speed k 
as a known information, we formulate the parameter estimation as an optimization problem. 

Problem 3.1. Assume that the occupation times Eij for i = 1,2, ... ,n and j = 1,2, ... ,m 
corresponding to a fixed time horizon T, different areas [a\, hW x [a^, 62] and various belt speeds 
Kj are given. Determine the parameters X,ai,a2 such that the least-squares deviation 



R{\,ai,a2) 



i=i j=i 



.p(A, cTi, CT2) - Eijj 



(5) 



is minimal. 



To minimize i?(A, di, (72) we use the simplex search method implemented in Matlab as 
the function fminsearch, see |5j . To test and demonstrate the parameter estimation using the 
optimization Problem |3.1| we proceed as follows: 



• Consider a process lA,cri,(T2 given values of the parameters A, ai and (T2- Compute the 
expected occupation times of the process corresponding to different areas [a\, b\] x [a^, b^] 
and belt speeds kj. 

• Check, if we can recover the given parameters A, a and a by solving the above Problem 

Let Eij be the expected occupation time of the process Yx=i^ai=i,a2=i corresponding 
to the time horizon T = 10, areas Vi = [0,1.0] x [-3.5,3.5], V2 = [0.5,1.5] x [-2.5,2.5], 
P3 = [1.0,2.5] X [-2.0,2.0] and P4 = [1.5,3.5] x [-1.25,1.25] and belt speeds ki = 1, k2 = 
2. We compute the expected occupation times Eij using a Monte-Carlo method based on 
5000 sample paths, see Table [2j Plugging this data into Problem 3.1 we obtain the cost 







V2 




P4 


Kl = 1 


1.23698 


1.16451 


1.63478 


1.92876 


K2 = 2 


0.81322 


0.70939 


0.93958 


1.06620 



Table 2: Expected occupation times. 

functional R{\, ai, a2). Solving the optimization problem we recover the following values for 
the parameters: 

A* = 0.98351, al = 1.01259, al = 1.00878 , 
which are quite close to the originnal set A = (Ti = (T2 = 1- 

Table [3] lists for some different settings the results of the numerical parameter estimation. 

Remark 3.1. Let < ai < 61 and < 02 < 62- Then, for a fixed belt speed k and time 
horizon T we have 



-^[ai,fei]x[a2,b2],K,T('^5Cri>'^2) — E[ai,bi] 



-a2_ 



which means, the expected occupation time distribution is symmetric around the central axis 
of the fiber lay-down, see Figure Due to this symmetry, we have just chosen symmetric 
domains T>i in the above examples. However, the selection of the areas should be adapted to the 
available fiber sample paths, since the occupation time distribution depends on the parameters 
of the lay-down process. 



true parameters 


recovered parameters 


A 


0-1 


-72 


A* 




0^2 


0.50 


1.00 


1.50 


0.50801 


1.00003 


1.45489 


0.60 


0.90 


1.40 


0.57202 


0.91942 


1.43237 


0.70 


2.00 


1.75 


0.68966 


2.10057 


1.77194 


1.00 


1.50 


1.30 


1.01343 


1.48292 


1.36154 


1.25 


1.25 


2.50 


1.23864 


1.28534 


2.52117 


1.40 


2.25 


0.80 


1.37193 


2.23118 


0.77826 


1.50 


0.80 


1.75 


1.53974 


0.84072 


1.73659 


2.00 


2.50 


2.50 


2.12084 


2.52146 


2.57379 


2.25 


2.50 


1.50 


2.11007 


2.47998 


1.51976 


2.50 


3.00 


3.00 


2.59092 


3.19728 


2.98174 



Table 3: List of recovered parameters. 



The recovered parameters in Table [s] have a maximal error of about 6%. This is a sufficient 
accuracy for many industrial applications. Nevertheless, increasing the number of sample 
paths used in the Monte-Carlo computations, we can improve the accuracy of the recovered 
parameters. 

For further testing the feasibility of our estimation method we use the fiber paths simulated 
by FYDIST, instead of the SDE @. On the left hand side of Figure [2| a FYDIST fiber path 
is shown. Based on such fiber paths we compute the expected occupation time for different 
areas Vi = [0.05,0.15] x [-0.039,0.039], V2 = [0.05,0.25] x [-0.025,0.025], V3 = [0.1,0.20] x 
[-0.018,0.018], P4 = [0.20,0.25] x [-0.01,0.01] and V5 = [0.25,0.29] x [-0.02,0.02], see 
Table [4} The belt speed ki = 0.0283 and the time horizon [0, 15.93] are fixed. 







V2 








K = 0.0283 


5.08290 


7.27732 


3.19642 


1.13889 


1.31337 



Table 4: Expected occupation times. 



Solving the optimization problem, we obtain the parameters 

A* = 4.328537, a* = 0.069968, = 0.029203 . (6) 

Figure [6] (top) shows a FYDIST fiber path used to calculate the Ei data and a sample path 
(bottom) of the fiber lay-down process iA*,CT*,(T* characterized by the estimated parameters 
and the known belt speed k = 0.0283. To compare the occupation time distribution of these 
two fiber paths, we consider the area [0, 0.46] x [—0.045, 0.045] and subdivide it into 20 vertical 
strips of uniform width, see Figure [6j Next, we compute the occupation time of both fiber 
paths with respect to each strip, see Figure [7} 
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Figure 6: Fiber path generated by FYDIST (top) and sample path of 1a*,ct*,(t| (bottom). 
2 p 1- -J* — I FYDIST^ — e — Process Y — ^ — Expected occupation time^-, — ^ ^-^ 
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Figure 7: Occupation times for the different strips and the two fiber paths shown in Figure 

m 

Conclusion 

Based on a hnear Ornstein-Uhlenbeck model ^ for the industrial fiber lay-down process, we 
have investigated the parameter estimation problem. For this model, we derived an analytical 
expression for the expected occupation time. This occupation time ([3| can be regarded as the 
mathematical equivalent to the mass per area density of the fiber web. Given the occupation 
time for different domains and different belt speeds, we were able to identify the parameters of 
the underlying Ornstein-Uhlenbeck process by solving the related least-squares minimization 
problem ([s]). Numerical computations based on Monte-Carlo simulations show the applicabil- 
ity of our method. We tested our method for fiber webs generated by an Ornstein-Uhlenbeck 
process as well as for webs generated using the industrial software FYDIST. Future research 
may focus on webs formed by more than one fiber and on nonlinear deposition models. 



9 



References 

[1] W. Bock, T. Gotz, M. Grothaus, and U. Liyanage. Parameter estimation from occupation 
times, 2011. arXiv:1105.5552vl [math.NA]. 

[2] L. Bonilla, T. Gotz, A. Klar, N. Marheineke, and R. Wegener. Hydrodynamic limit of 
a Fokker-Planck equation describing fiber lay-down processes. SI AM J. Appl. Math., 
68(3):648-665, 2008. 

[3] T. Gotz, A. Klar, N. Marheineke, and R. Wegener. A stochastic model for the fiber 
lay-down process in the nonwoven production. SIAM J. Appl. Math., 67(6):1704-1717, 
2007. 

[4] P. Kloeden and E. Platen. Numerical Solution of Stochastic Differential Equations. 
Springer-Verlag, 1992. 

[5] MathWorks. MATLAB 7 Function Reference, 2007. 

[6] T. Miiller-Gronbach and K. Ritter. Variable subspace sampling and multi-level algorithms. 

In P. L'Ecuyer and A. Owen, editors, Monte Carlo and Quasi-Monte Carlo Methods 2008, 
pages 131-156. Springer-Verlag, Berlin, 2009. 



10 



